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ABSTRACT 

Aims. This work presents a uniform and homogeneous study of chemical abundances of refractory elements in 101 stars with 
and 94 without known planetary companions. We carry out an in-depth investigation of the abundances of Si, Ca, Sc, Ti, V, 
Cr, Mn, Co, Ni, Na, Mg and Al. The new comparison sample, spanning the metallicity range —0.70 < [Fe/H] < 0.50, fills the 
gap that previously existed, mainly at high metallicities, in the number of stars without known planets. 

Methods. We used an enlarged set of data including new observations, especially for the field “single” comparison stars . The 
line list previously studied by other authors was improved: on average we analysed 90 spectral lines in every spectrum and 
carefully measured more than 16 600 equivalent widths (EW) to calculate the abundances. 

Results. We investigate possible differences between the chemical abundances of the two groups of stars, both with and without 
planets. The results are globally comparable to those obtained by other authors, and in most cases the abundance trends of 
planet-host stars are very similar to those of the comparison sample. 

Conclusions. This work represents a step towards the comprehension of recently discovered planetary systems. These results could 
also be useful for verifying galactic models at high metallicities and consequently improve our knowledge of stellar nucleosynthesis 
and galactic chemical evolution. 
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1. Introduction 

Over the last ten years a large number of stars harbouring 
planets have been found. The first giant planet was discov¬ 
ered around 51 Peg by Mayor & Queloz (1995) and there 
are more than 150 planetary-nrass companions presently 
known orbiting solar-type stars. The growing number of 
extrasolar planets 1 has activated intensive study of these 
objects and their parent stars, and nowadays extensive 


Send offprint requests to: G. Israelian, e-mail: gil@iac.es 

* Based on observations collected at the La Silla Observatory, 
ESO (Chile), with CORALIE spectrograph at the 1.2 m Euler 
Swiss telescope, and with the FEROS spectrograph at the 1.52 
m ESO telescope, at the Paranal Observatory, ESO (Chile), 
using the UVES spectrograph at the VLT/UT2 Kueyen tele¬ 
scope, and with the UES and SARG spectrographs at the 4 
m William Hershel Telescope (WHT) and the 3.5 m TNG, re¬ 
spectively, both at La Palma (Canary Islands). 

1 A complete updated table of known planets can be found 
at https://obswww.unige.ch/exoplanets. 


studies of the properties of the new planetary systems are 
conceivable. Several spectroscopic analyses of iron abun¬ 
dances (Gonzalez 1998; Gonzalez et al. 2001; Laws et 
al. 2003; Santos et al. 2001a,2001b, 2003, 2004a,b, 2005) 
have suggested that planet-host stars are more metal-rich 
than field dwarfs. These results show that the probabil¬ 
ity of finding a planet is a strongly increasing function of 
stellar metallicity, at least for [Fe/H] above solar value. 
Two interpretations, the self-enrichment hypothesis and 
the primordial hypothesis , have been proposed to explain 
a possible connection between the metallicity excess and 
the presence of planets. The former considers that the 
observed iron overabundances derives from the accretion 
of a large amount of rocky planetesimal material on to 
the star (Gonzalez 1997). The latter, the “primordial sce¬ 
nario” , suggests that the iron excess in stars with planets 
just reflects the high metal content of the protoplanetary 
cloud from which stars and planets were formed (Santos 
et al. 2000, 2001a,2001b). 
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In this context, abundance trends of chemical species 
other than iron can give important clues in this debate, 
so that discriminating between these two possibilities will 
help in understanding how planetary systems form. Efforts 
have been made to analyse the chemical abundances of 
light elements (e.g. Garcia Lopez & Perez de Taoro 1998; 
Gonzalez & Laws 2000; Ryan 2000; Deliyannis et al. 2000; 
Israelian et al. 2003a, 2004; Santos et al. 2002, 2004b), as 
well as other metals (e.g. Gonzalez & Laws 2000; Gonzalez 
et al. 2001; Smith et al. 2001; Takeda et al. 2001; Sadakane 
et al. 2002; Fisher & Valenti 2005). Most of these au¬ 
thors considered inhomogeneous comparison samples of 
field dwarfs from the literature that might be a source 
of systematic errors. The majority of these studies sup¬ 
port the primordial hypothesis (Pinsonneault et al. 2001; 
Santos et al. 2001a,2001b, 2002; Sadakane et al. 2002), 
but evidence of pollution has been found for a few cases 
(Israelian et al. 2001,2003a,2003b; Low & Gonzalez 2001). 

Unbiased and homogeneous studies of Fe abundance 
in stars with and without planets have been performed 
by Santos et al. (2001a, 2003, 2004a, 2005). Similar stud¬ 
ies have been recently carried out for elements other than 
iron. Refractory elements (some a and Fe-group elements) 
have been analysed by Bodaghee et al. (2003); all volatile 
elements (C, S, Zn, N, O) have been studied by Ecuvillon 
et al. (2004a,2004b, 2005) and Na, Al and Mg by Beirao 
et al. (2005). These studies all required a uniform high- 
metallicity comparison, given the lack of “single” field 
stars in the data for [Fe/H]> 0.1. We provide results us¬ 
ing new and more precise atmospheric parameters from 
high quality spectra and also complete the high-metallicity 
comparison between the two samples. 

In this paper we present a detailed, homogeneous and 
uniform study of Si, Ca, Ti, Sc, V, Cr, Mn, Co, Ni, Na, 
Mg and Al in a set of 101 planet-harbouring stars and a 
group of 94 solar-type stars with no known planets in the 
metallicity range —0.70 <[Fe/H] < 0.45. We have made 
use of about 16 600 EWs measured to calculate the [A'/H] 
ratio for each element and have plotted the results in the 
[X/H] vs. [Fe/H] and [X/Fe] vs. [Fe/H] planes. 

2. Data 

All the objects from the comparison sample belong to the 
CORALIE extrasolar-planets-finding programme 2 . The 
high resolution spectra are the same as those used by 
Santos et al. (2004a) to derive precise and uniform stel¬ 
lar parameters for 98 planet-host stars and 41 compar¬ 
ison sample “single” dwarfs in a volume-limited sample 
in the solar neigborhood (< 20 pc). It should be pointed 
out that the star HD 219542 B was excluded from the 
planet-host list presented in Santos et al. (2004a) since 
the presence of a planet around this star has been rejected 
(Desidera et al. 2004). All our spectra were collected dur¬ 
ing several observing campaigns with the CORALIE spec¬ 
trograph on the 1.2 nr Euler Swiss Telescope, the FEROS 

2 http://unige.ch/ udry/planet/planet.html. 


Table 1 . Atomic parameters of the spectral lines used for 
every element. Col. 1: wavelength (in A); Col. 2: excitation 
energy of the lower energy level in the transition (in eV); Col. 3: 
oscillator strengths based on an inverse solar analysis. 


A 

Xi 

log gf 

A 

Xl 

log gf 

Si I ; log e 

O = 7.55 A = 14 

Cr I ; log < 

= 5.67 A = 24 

5665.56 

4.92 

-1.980 

5304.18 

3.46 

-0.680 

5690.43 

4.93 

-1.790 

5312.86 

3.45 

-0.580 

5701.10 

4.93 

-2.020 

5318.77 

3.44 

-0.710 

5772.14 

5.08 

-1.620 

5480.51 

3.50 

-0.830 

5793.09 

4.93 

-1.910 

5574.39 

4.45 

-0.480 

5948.55 

5.08 

-1.110 

5783.07 

3.32 

-0.400 

6125.02 

5.61 

-1.520 

5783.87 

3.32 

-0.150 

6142.49 

5.62 

-1.480 

5787.92 

3.32 

-0.110 

6145.02 

5.61 

-1.400 

Mn I ; log 

e Q = 5.39 A = 25 

6155.15 

5.62 

-0.750 

4265.92 

2.94 

- 0.440 

6721.86 

5.86 

-1.090 

4470.13 

2.94 

- 0.550 

Ca I ; log e Q = 

6.36 A = 20 

4502.13 

2.92 

- 0.490 

5512.98 

2.93 

-0.440 

5399.47 

3.85 

-0.0969 

5581.97 

2.52 

-0.650 

5413.68 

3.85 

-0.470 

5590.12 

2.52 

-0.710 

5432.54 

0.00 

-3.620 

5867.56 

2.93 

-1.590 

6440.93 

3.77 

-1.250 

6161.29 

2.52 

-1.220 

Co I ; log 

e 0 = 4.92 A = 27 

6166.44 

2.52 

-1.120 

5301.04 

1.71 

-1.930 

6169.05 

2.52 

-0.730 

5325.27 

4.02 

-0.120 

6169.56 

2.52 

-0.440 

5342.70 

4.02 

0.574 

6449.82 

2.52 

-0.630 

5483.36 

1.71 

-1.220 

6455.60 

2.52 

-1.370 

5647.23 

2.28 

-1.580 

Sc II ; log 

e D = 

3.10 A = 21 

6093.15 

1.74 

-2.340 

5239.82 

1.45 

-0.760 

6455.00 

3.63 

-0.280 

5318.36 

1.36 

-1.700 

Ni I ; log f 

: 0 = 6.25 A = 28 

5526.82 

1.77 

0.150 

5578.72 

1.68 

-2.650 

6245.62 

1.51 

-1.040 

5587.86 

1.93 

-2.380 

6300.69 

1.51 

-1.960 

5682.20 

4.10 

-0.390 

6320.84 

1.50 

-1.840 

5694.99 

4.09 

-0.600 

6604.60 

1.36 

-1.160 

5805.22 

4.17 

-0.580 

Ti I ; log c 

: Q = 4.99 A = 22 

5847.00 

1.68 

-3.410 

5471.20 

1.44 

-1.550 

6086.28 

4.26 

-0.440 

5474.23 

1.46 

-1.360 

6111.07 

4.09 

-0.800 

5490.15 

1.46 

-0.980 

6128.98 

1.68 

-3.370 

5866.46 

1.07 

-0.840 

6130.14 

4.26 

-0.950 

6091.18 

2.27 

-0.460 

Na I ; log 

£o = 6.33 A = 11 

6126.22 

1.07 

-0.410 

5688.22 

2.104 

- 0.625 

6258.11 

1.44 

-0.440 

6154.23 

2.102 

- 1.607 

6261.11 

1.43 

-0.490 

6160.75 

2.104 

- 1.316 

6303.76 

1.44 

-1.600 

Mg I ; log 

e 0 = 7.58 A = 12 

6312.24 

1.46 

-1.580 

5711.09 

4.346 

- 1.706 

V I ; log e, 

= = 4.00 A — 23 

6318.72 

5.108 

- 1.996 

5668.37 

1.08 

-1.00 

6319.24 

5.108 

- 2.179 

5670.85 

1.08 

-0.460 

8712.69 

5.932 

- 1.204 

5727.05 

1.08 

-0.000 

8736.02 

5.946 

- 0.224 

5727.66 

1.05 

-0.890 

Al I ; log e 

: 0 = 6.47 A = 13 

5737.07 

1.06 

-0.770 

6696.03 

3.143 

- 1.570 

6090.21 

1.08 

-0.150 

6698.67 

3.143 

- 1.879 

6216.35 

0.28 

-0.900 

7835.31 

4.022 

- 0.728 

6285.16 

0.28 

-1.650 

7836.13 

4.022 

- 0.559 

6452.31 

1.19 

-0.820 

8772.87 

8773.91 

4.022 

4.022 

- 0.425 

- 0.212 


spectrograph on the 2.2 m ESO/MPI Telescope (both at 
La Silla, Chile), the UVES spectrograph on the VLT/UT2 
Kueyen Telescope (Paranal Observatory, ESO, Chile), the 
SARG spectrograph on the 3.5 m TNG and the UES spec¬ 
trograph on the 4.2 m WHT (both at the Roque de los 
Muchachos Observatory, La Palma, Spain). The spectro¬ 
graph and data characterisctics are listed in Table |2j 
In order to extend the number of stars in a comparison 
sample, we added 48 FEROS spectra from Santos et al. 
(2005) and five new SARG spectra. We refer the reader 
to this paper for a description of the data. The S/N ratio 
of the spectra varies between 150 and 350. New high S/N 
spectra obtained in 2004 with VLT/UVES have been used 
as well (see Ecuvillon et al. 2004a,2004b, 2005) 
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Table 2. Spectrographs and data characteristics. 


Spectrograph/Telescope 

Observatory 

Resolution 

(A/AA) 

Range 

(A) 

CORALIE/1.2-m Ewler Swiss 

La Silla(Chile) 

50000 

3800-6800 

FEROS/1.52-m ESO 

La Silla (Chile) 

48000 

3600-9200 

UVES/VLT 8-mKuyen UT2 

Paranal (Chile) 

110000 

4800-6800 

SARG/3.5-m TNG 

ORM (la Palma) 

57000 

5100-10100 

UES/4-m WHT 

ORM (la Palma) 

55000 

4600-7800 



1500 5000 5500 6000 4500 5000 5500 6000 4500 5000 5500 6000 

T_eff T_eff T_eff 



3.8 4.0 4.2 4.4 4.6 4.8 5 .( 
I°g g 


4.0 4.2 4.4 4.6 4.8 3.8 4.0 4.2 4.4 4.6 4.8 5.0 
log g log g 


Fig. 1. [X /Fe] vs. T e g and [A'/Fe] vs. logg. Triangles and asterisks are planet host and comparison sample stars, 
respectively. Dotted lines represent the solar value. The slopes represent linear least-squares fits to both planet-hosts 
and comparison stars. 


3. Spectral analysis 

Chemical abundances of the nine refractory elements stud¬ 
ied here, as well as Na, Mg and Al, were derived perform¬ 
ing a standard local thermodynamic equilibrium (LTE) 
analysis, strictly differential with respect to the Sun. Solar 
abundances for each element were taken from Anders & 
Grevesse (1989) using a solar model with T e g = 5777 K, 
log<? = 4.44 dex, = 1.0 km s . We used a revised ver¬ 
sion (2002) of the MOOG code (Sneden 1973) 3 (with the 
abfind driver) and a grid of Kurucz (1993) ATLAS9 at¬ 
mospheres. The atmospheric parameters, effective temper¬ 
ature (T e ff), surface gravity (log g), microturbolence (£ t ) 
and metallicity ([Fe/H]), and their corresponding uncer¬ 
tainties, were taken from Santos et al. (2004a, 2005). It 
should be stressed that these stellar parameters were de¬ 
rived in a uniform way, using high resolution spectra, the 
same line list and model atmospheres for all the stars. The 
spectral lines of the refractory elements analysed here were 
extracted from the study by Bodaghee et al. (2003), while 
the line lists of Na, Mg and Al were taken from Beirao 
et al. (2005). These lists were successively modified in or¬ 
der to minimize abundance errors. Since the V I line at 
6531.42 A, the Mn I line at 5388.50 A and the Co I lines at 

3 The MOOG2002 source code can be downloaded at 
http://verdi.as.utexas.edu/moog.html 


5312.86 A and 6632.44 A were difficult to measure in most 
of the spectra (too weak or blended), we eliminated them 
from the list. We also excluded the Mg line at 8923.57 A 
because it was not measurable in the new Feros spectra 
analysed for Na, Mg and Al. We added more lines of V, 
Mn and Co from Gurtovenko & Kostik (1989). Before in¬ 
cluding these lines in our list, we first verified that each 
line was not too strong and checked for possible blend¬ 
ing, using the Kurucz Solar Atlas (Kurucz et al. 1984). 
For these new lines we derived semi-empirical atomic os¬ 
cillator strengths using their EWs measured in the solar 
atmosphere with (T e g, log g, £*)= 5770 K, 4.44 dex and 1.0 
km s -1 ), and performed an inverted solar analysis. The 
number of lines increased from 5 to 9 for V, from 3 to 8 
for Mn, and from 5 to 7 for Co. 

Finally, we considered about 80 spectral lines for the 
analysis of refractory elements in each spectrum and about 
14 lines for the study of Na, Mg and Al abundances (see 
Table nj in 53 new FEROS and SARG spectra. We have 
measured about 16 600 equivalent widths in our spectra. 
For each spectral line EWs were measured by a Gaussian 
fit using the splot task within the “echelle” IRAF pack¬ 
age. 4 Furthermore we employed a new program (IRAF 

4 IRAF is distributed by National Optical Astronomy 
Observatories, operated by the Association of Universities for 

























Relative Freq 
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Distribution T_eff 



Fig. 2. Distributions of planet-host stars and non-planet- 
host stars as functions of T e g. The planet-host stars and 
the comparison sample are represented with solid and dot¬ 
ted lines, respectively. 


Distribution [Ti/H] Distribution [V/H] Distribution [Co/H] 



Fig. 3. [X/H] distributions for Ti, V and Co for the two 
subsamples of planet host (solid line) and comparison sam¬ 
ple (dotted lines) stars with the same T e g distributions. 


scripts) written by N. C. Santos for automatic measure 
of the EWs of Na, Mg and Al lines. This program was 
previously used to measure the EWs of a list of Fe lines 
(Santos et al. 2005). Every measurement was taken care¬ 
fully achieving the best agreement between lines profile 
and Gaussian fits. About 25% of the equivalent widths 
used to calculate abundances come from Bodaghee et 
al. (2003), whose stellar parameters have been updated by 
Santos et al. (2004a, 2005). Detailed observational data 
(e.g. line-by-line EWs for each element) are available as 
an electronic table (Table 11) at the Centre des Donnees 
Strasbourg (CDS). 


Research in Astronomy, Inc., under contract with the National 
Science Fundation, USA. 


Table 3. Slopes of [X /Fe] ratios as functions of effective 
temperatures per 1000 K (in dex) for both samples (planet 
hosts and non-hosts). 


Species 

Slopes ± rms 

Species 

Slopes ± rms 

Si 

-0.015 

± 

0.010 

Cr 

-0.007 

± 

0.007 

Ca 

0.057 

± 

0.013 

Mn 

-0.029 

± 

0.020 

Sc 

-0.014 

± 

0.011 

Co 

-0.122 

± 

0.012 

Ti 

-0.137 

± 

0.015 

Ni 

0.023 

± 

0.008 

V 

-0.242 

± 

0.015 






3.1. Uncertainties 

We have tested the dependence of our results on atmo¬ 
spheric parameters (Fig. |TJ) and refer the reader to the 
Figs 1 -2 by Beirao et al. (2005) for the plots of [. X /Fe] 
( X = Na, Mg, Al) vs. T e g and vs. log g. In Table 0 are 
listed the slopes of the [X/Fe] ratios for refractory ele¬ 
ments as a function of T e g for all stars in the distribu¬ 
tions. We observe no characteristic trends for the major¬ 
ity of the elements studied. In the case of Ti, V and Co 
we note decreasing trends for the [X/Fe] ratio correspond¬ 
ing to T e ff values greater than 5500 K. This trend could be 
caused by NLTE effects which, however, contribute mostly 
in the case of the more metal-poor stars. We do not take 
NLTE effects into account because they are usually small 
in metal-rich stars for complex atoms such as iron (see 
Edvardsson et al. 1993; Thevenin & Idiart 1999). On the 
other hand, although NLTE effects have not been stud¬ 
ied for a number of elements (e.g. V, Co, etc.), they are 
nonetheless not thought to alter the main conclusion of 
the present paper. Some unknown blended lines may also 
be responsible for these trends. In fact, the overestimation 
of EWs caused by the increasing blending effects becomes 
more severe at lower T e g. 

The distribution of planet-hosts and non-planet-host 
stars as function of T e g shows that the latter are on aver¬ 
age cooler (see Fig. [21 . To verify this, we selected randomly 
two subgroups of stars (planet and no-planet hosts) with 
the same T e g distributions and obtained the [X/H] distri¬ 
bution for Ti, V and Co. We repeated this process for 200 
randomly selected subgroups and calculated the average 
distribution for each element (Ti, V, Co). In Fig.Olwe ob¬ 
served the same behaviour as that found for all the stars 
analysed. In both cases the differences between average 
[X'/H] for hosts and no-hosts are of the same order (see 
Table [21 

In contrast, systematic errors are difficult to locate but 
are largely reduced by good data quality and good instru¬ 
mental resolution. Since we have analysed more than one 
line for each element (usually 7 8 lines) in a given star, 
the total dispersions around the average abundance are 
more significant compared to the continuum observed un¬ 
certainties, which are usually around 0.05 dex. 

When analysing many lines to calculate abundances, 
uncertainties in the atmospheric parameters should be the 
primary source of abundance errors. These are of the or- 
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Table 4. Sensitivity of the abundances of Si, Ca, Sc, Ti, V, Cr, Mn, Co and Ni to changes of 100 K in T e g, 0.30 dex in logg 
and [Fe/H], 0.50 km s -1 in £t. 



der of 50 K in T e g, 0.12 dex in logg, 0.08 km s -1 in 
and 0.05 dex in [Fe/H] (see Santos et al. 2004b, 2005). 
The abundance sensitivity to changes in atmospheric pa¬ 
rameters were estimated as follows. First we selected a 
set of 12 stars from our list: for each atmospheric pa¬ 
rameter we chose three stars with similar values for all 
the parameters except for the one considered, which must 
vary within the sample. We then generated new atmo¬ 
spheric models, changing only the “varying” parameter 
and calculating the corresponding abundances. The dif¬ 
ference between these values and those obtained without 
varying the parameter, give us the abundance sensitiv¬ 
ity to changes in the parameters. Tables 0 and El show 
these results after varying the effective temperature by 
A T e g = ±100 , surface gravity by Alogg = ±0.30 dex, 
metallicity by A [Fe/H] = ±0.30 dex and microturbulence 
by A£* = ±0.50ferns -1 . We note that ions such as Sc II 
are generally more sensitive to changes in surface gravity, 
while neutral atoms are influenced mostly by uncertain¬ 
ties in the effective temperature. At a glance, we observe 
from TableEJthat the V and Ti abundances can vary with 
temperature changes of A[X/H] ~ ±0.12 dex, while we 
associate A [Sc/H] ~ ±0.11 dex with changes in surface 
gravity. With respect to Table 0 we note that Na abun¬ 


dances are more sensitive to variations in atmospheric pa¬ 
rameters (e.g. T e ff and log g) than Mg and Al. Finally, 
we evaluated the errors in the abundances of all the el¬ 
ements, adding quadratically the standard deviation of 
the mean abundance obtained from all the measured lines 
and the uncertainties dued to the abundance sensitivities 
to changes in the atmospheric parameters. For each star 
these ’’total” errors are of the order of 0.10 dex. 

4. Results 

Bodaghee et al. (2003) recently carried out a spectroscopic 
analysis of the same refractory elements as those pre¬ 
sented in our work while Beirao et al. (2005) did the same 
for Na, Mg and Al. These authors did not find any sig¬ 
nificant differences between planet host and comparison- 
sample stars, a result in perfect agreement with our find¬ 
ings. Because of the lack of comparison sample stars with 
[Fe/H] >0.1 dex in previous studies, the abundance dis¬ 
tribution of stars with giant planets looked like a high 
metallicity extension of the curves traced by field dwarfs 
without planets. New spectra of metal-rich stars with no 
planets have been gathered, and consequently a complete 
comparison is also possible in the high metallicity domain. 
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Table 5. The same as Table^Jbut for Na, Mg and Al. 



It is still conceivable that certain specific trends found for 
the metal-rich tail of Galactic chemical evolution are due 
to the presence of planets (see plots for Co, V, Na, Mg, Al 
in Figs mm However, we should stress that even the new 
comparison sample contains about 20% of stars at [Fe/H] 
>0.1. All the results, together with their total errors (see 
Sec. Id.ID . are listed in the Tables 7-10. Furthermore, we 
have compared the abundances calculated with EW values 
measured with different instruments. As shown in Fig.QJ 
we did not detect any significant discrepancies among the 
abundance values for the nine refractory elements. 

The histograms provide the distributions of [AT/H] 
with A = Si, Ca, Sc, Ti, V, Cr, Mn, Co and Ni (Fig.EJ), 
and with X = Na, Mg, Al (Fig. for the two samples 
of stars, with and without planets. These results are sim¬ 
ilar to those presented for [Fe/H] by Santos et al. (2004a, 
2005) and clearly confirm that the observed metallicity 
excess is, as expected when extended to elements other 
than iron. We observe similar features to those noted in 
the [A/H] distribution for refractory elements (Bodaghee 
et al. 2003) and for iron (Santos et al. 2001a, 2004a; Reid 
2002). For example, these histograms of planet-host stars 
are usually not symmetrical. This interesting feature is 
particularly evident for Ca, Sc, Co, Ni, Na and Al, for 
which the distribution seems to be an increasing func¬ 


tion of [A'/H] up to a certain value, after which it falls 
abruptly. This cut-off corresponds to [A'/H] ~ 0.5 for Si, 
Ti, Ni, Na, [A/H] ~ 0.7 for V, Cr, Mn and Co and [A'/H] 
~ 0.6 for Mg and Al. Only in the case of Ca does the 
distribution fall to [A/H] ~ 0.3. For some elements (e.g. 
Ti, V, Cr, Mg) the distributions appear to be slightly bi- 
modal. This is probably related to a lack of stars with 
[A/H] ~ 0.3 in the planet-host sample and [A/H] ~ —0.2 
in the comparison group for these elements. The average 
values <[A/H]>, the rms dispersions for the two distri¬ 
bution, and the difference between the average [A'/H] for 
stars with and without planets are listed in Table EH We 
note that this difference varies from 0.13 dex for V to 0.29 
dex for Mn and Na. The difference between the average 
abundance values of the two groups (see Table 0 is only 
an estimate and is not very significant, given the usually 
high dispersion around the mean values. 

In Fig. □ are shown [A/H] vs. [Fe/H] plots for re¬ 
fractory elements. The [A/H] ratio is a linear function 
of [Fe/H] and the small amount of scatter give a certain 
plausibility to our results. We refer the reader to the Fig. 5 
by Beirao et al. (2005) for the plots of [A'/H] vs. [Fe/H] 
(X= Na, Mg and Al). These plots are not presented here 
because adding 53 new comparison sample stars did not 
change the results for these elements. However, [A/Fe] vs. 
[Fe/H] plots for Na, Mg and Al are shown in Figure ITol 

5. Galactic chemical evolution trends 

With the exception of the lightest elements (e.g. H and 
He), the history of the Galaxy’s chemical composition is 
dominated by nucleosynthesis occurring in many genera¬ 
tions of stars (McWilliam 1997). The low-mass stars are 
like “fossils” because their lifetimes are sometimes compa¬ 
rable to the age of the Galaxy. It might actually be sup¬ 
posed that, at least for F-G dwarfs, the external envelope 
of the stars have preserved much of their original chemi- 

Table 6. Average abundances <[A'/H]> for stars with planets 
and for comparison sample stars. For each element the rms 
around the mean and the abundances difference between the 
two samples (101 and 94 respectively) are also listed. 


Species 

(X) 

Planet-hosts 
<[X/H]> a 

No-hosts 
<[X/H]> a 

Average 

diff. 

Si 

0.09 

0.20 

-0.13 

0.24 

0.23 

Ca 

-0.01 

0.19 

-0.22 

0.22 

0.22 

Ti 

0.15 

0.17 

0.01 

0.19 

0.14 

Sc 

0.16 

0.20 

-0.08 

0.24 

0.25 

V 

0.21 

0.24 

0.08 

0.26 

0.13 

Cr 

0.04 

0.22 

-0.18 

0.25 

0.22 

Mn 

0.15 

0.33 

-0.14 

0.33 

0.29 

Co 

0.16 

0.25 

-0.04 

0.26 

0.20 

Ni 

0.09 

0.24 

-0.17 

0.27 

0.26 

Na 

0.12 

0.24 

-0.17 

0.26 

0.29 

Mg 

0.19 

0.19 

-0.05 

0.21 

0.23 

Al 

0.19 

0.21 

-0.01 

0.23 

0.18 
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Fig. 4. Comparison of the abundances of nine refractory elements calculated with EW values measured from spectra taken 
with different instruments. Along the rr-asis are results from FEROS spectra, while along the y-axis we represent abundance 
values for the same objects, calculated using other spectra. UVES, CORALIE and SARG results are respectively represented 
with rhornbi (HD 16141, HD 27442), triangles (HD 83443, HD 92788, HD 106252, HD 114386) and asterisks (HD 114783, 
HD 147513). 








Fig. 5. [X /H] vs. [Fe/H] for Si, Ca, Ti (left), for Sc, Mn and V (centre), for Cr, Co and Ni (right). Triangles and 
asterisks represent stars with planetary-mass companion and “single” field stars, respectively. The intersection of the 
dotted lines indicates solar value. 
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Fig. 6. [X /H] distribution for each element. The planet-host stars and the comparison sample are represented with 
solid and dotted lines, respectively. 


cal composition, since it has not been convectively mixed 
with internal matter. The relation [A/Fe] vs. [Fe/H] is 
traditionally used in observational studies of the chemical 
evolution of the Galaxy since iron is a good chronological 
indicator of nucleosynthesis (controversial). Furthermore, 
iron lines are numerous and easy to measure in the spectra 
of dwarfs. 

In Fig. m we present [A'/Fe] vs. [Fe/H] plots for Na, 
Mg and Al taken from the work by Beirao et al. (2005) 
with the addition of abundances from a new comparison 
sample. We also did calculations relative to the iron abun¬ 
dances ([A'/Fe]) of a elements (e.g. Si, Ca, Ti) and of 
iron-peak elements (e.g. Cr, Mn, Co, Ni). The former are 
believed to be mostly produced in the aftermath of ex¬ 
plosions of type II supernovae (SNe II), although, follow¬ 
ing some models, these elements might also be produced 
during a type la supernova (SN la) event (Thielmann et 
al. 2002). Meanwhile, most of the latter would have been 
synthesized by SNe la explosions. Magnesium is supposed 
to be produced by SNe II, thus comparing Mg abundances 
with those found for other a elements could provide us 
with evidence concerning the origin of these elements. 
Sodium and aluminium are thought to be mostly a prod¬ 


uct of Ne and C burning in massive stars. It is not clear 
how Sc is formed, because in the periodic table it is in¬ 
termediate between a and iron-peak elements. The origin 
of manganese is also debated. However, Sc and Mn abun¬ 
dances in long-lived F and G stars are of great interest 
since they could also introduce special constraints on nu¬ 
cleosynthesis theory (Nissen et al. 2000). Even though the 
main aim of this study was to compare the abundances of 
refractory elements in stars with and without planets, our 
results also give us a chance to increase our present knowl¬ 
edge of the chemical evolution of the Galaxy at [Fe/H] > 0. 
In FigsQ and m we present the [A/Fe] vs. [Fe/H] plots 
for all the elements. 

6. Comparison with the literature 

There are already several studies concerning the chemical 
abundances of elements other than iron in F, G, K main 
sequence stars in the solar neighborhood (see Edvardson 
et al. 1993; Feltzing & Gustafsson 1998; Chen et al. 2000; 
Nissen et al. 2000; Fulbright et al. 2002; Reddy et al. 2003; 
Bensby et al. 2003; Allende Prieto et al. 2004). In addi¬ 
tion, various studies on chemical abundances in planet- 
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Fig. 7. [X/Fe] vs. [Fe/H ](in dex ) for a elements (Si, Ca, Ti, Sc) and for iron-peak elements (Cr, Co, Ni, Mn, V). 
[V/Fe] vs. [Fe/H] for T e s > 5000 K objects only (below right panel). The triangles and asterisks represent stars with 
planets and “single” stars, respectively. The intersection of the dotted lines indicates the solar value. 


host stars have gradually emerged (Gonzalez et al. 2000; et al. 2005; Fischer & Valenti 2005). We note that our 
Santos et al. 2000; Takeda et al. 2001; Sadakane et al. 2002; abundance trends generally agree with those published in 
Bodaghee et al. 2003; Ecuvillon et al. 2004a, 2004b; Beirao the literature and this lends a certain reliability to our re- 
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suits. In the following subsections, we describe [ X /Fe] vs. 
[Fe/H] trends for each element and then make a brief com¬ 
parison with studies in the literature on this subject. This 
comparison is divided into two parts: first, we compare 
studies on abundances in stars with planets; second, we ex¬ 
tend the comparison to other studies regarding abundance 
trends in metal-rich stars of the Galactic disc. We focus 
our attention mainly on the [Fe/H] > 0 range. Given the 
unobserved or probably insignificant differences between 
the two samples of stars presented, we consider the distri¬ 
bution as a whole for the rest of the analysis. Possible dif¬ 
ferences between our trends and those recently published 
could be of great interest since it might reflect the pres¬ 
ence of planets. In contrast, any relation to the presence 
of a planet is perhaps coincidental and the global trends 
observed are probably best interpreted as a consequence 
of Galactic chemical evolution. 

6.1. The a elements 

6.1.1. Silicon 

In Figs 0 and 0 we note a slight [Si/Fe] overabundance 
compared to the solar value at [Fe/H] < —0.2 (with a 
rather larger scatter), while the point distribution remains 
constant at [Si/Fe] ~ 0.0 for the rest of the metallicity 
range. 

Our results agree with the previous chemical analyses 
of stars with planets by Gonzalez et al. (2000, hereafter 
GZ00), Bodaghee et al. (2003, hereafter BOD) and Fischer 
& Valenti (2005). Sadakane et al. (2002, hereafter SD02) 
in their study analysed only a few objects (12 stars with 
extrasolar planets) but their abundance distribution is flat 
at high metallicities. 

Edvardsson et al. (1993, hereafter EAGLNT) and 
Chen et al. (2000, hereafter COO) obtained similar re¬ 
sults, as well as Bensby at al. (2003, hereafter BEN) and 
Fulbright et al. (2002, hereafter FUL) in the range [Fe/H] 
> 0. In another study on metal-rich stars ([Fe/H] > 0.10) 
(Felzing & Gustafsson 1998, hereafter FG98) the [Si/Fe] 
distribution exhibits a constant trend around the solar 
value. All these results differ from those of Allende Prieto 
et al. (2004, hereafter AL04), in which the abundance 
trend for Si changes abruptly around [Fe/H] ~ 0 and as¬ 
sumes a positive slope. 

6.1.2. Calcium 

Contrary to the other a elements, the [Ca/Fe] ratio seems 
to decrease quite uniformly (see Fig.0and Fig. 0 ). In par¬ 
ticular, the plot suggests the presence of a plateau in the 
range —0.2 < [Fe/H] < 0.2 followed, for higher metallici¬ 
ties, by a slight fall-off. Planet-host stars with T e ff < 5000 
K show a certain dispersion in [Ca/Fe] values, for example 
the stars HD 177830, HD 137759 and HD 114783. 

The similar behaviour of calcium trends for stars both 
with and without planets has been observed in other stud¬ 
ies of chemical abundances (BOD and GZ00). As BOD 


have noticed, SD02 analysed only two stars in the range 
0.2 < [Fe/H] < 0.4 (where a possible plateau is suggested), 
so the comparison with our results is only partially valid. 

The Ca distribution appears quite “flat” in the data 
presented by EAGLNT, FG89, AL04 and FUL. However, 
the results plotted by BEN appear to decrease for [Fe/H] 
~ 0.2 -0.3. 

6.1.3. Titanium 

The titanium [Ti/Fe] ratio decreases in the range —0.6 < 
[Fe/H] < 0 until the solar value where the distribution 
settles (Fig. 0 and 0 . 

Similar trends have been obtained by BOD and Fischer 
& Valenti (2005), while we cannot observe any slope 
change around —0.2 < [Fe/H] < 0 in SD02, because only 
a few objects have been plotted. Despite the large scatter 
in the GZ00 data, titanium abundances seem to decrease 
gradually with metallicity. 

A quite pronounced point dispersion is also observed in 
EAGLNT, BEN and CH00, especially for [Fe/H] < -0.3 
(as seen here). This scatter is probably due to the overes¬ 
timation of EWs caused by the increasing blending effect 
(which becomes severer as T e g falls) or by Galactic evolu¬ 
tion effects. Finally, comparing the trends for [Fe/H] > 0 
we note that [Ti/Fe] values remain approximately con¬ 
stant in the plots of BEN, FUL and FG98, while the data 
of AL04 again show a rise above solar metalicity. 

6.1.4. Scandium 

Figures 0 and 0 illustrate that the scandium trend is sim¬ 
ilar to that of other a elements. It drops until —0.2 < 
[Fe/H] < 0 and afterwards the [Sc/Fe] ratio approximately 
follows the solar value at high metallicities. 

Previous studies on stars with planets have obtained 
similar results, as presented in BOD and SD02. 

There are only a few studies in the literature about 
scandium abundances. The most detailed one was pre¬ 
sented by Nissen et al. (2000) in the range —1.4 < [Fe/H] 
< 0.1, so comparison with our results at high metallicity 
is not possible at all. However, it is interesting to note a 
slope change at [Fe/H] ~ —0.3, —0.2 that is also repre¬ 
sented in our graphs. Results by FG98 show a substantial 
star-to-star scatter but globally the [Sc/Fe] ratio remains 
around solar for [Fe/H] > 0. In contrast, the AL04 plot 
shows large enrichments in scandium for metal-rich stars. 

6.2. The Fe-group elements 
6.2.1. Manganese 

Manganese is one of the lesser studied element in the lit¬ 
erature. Manganese lines turn out to be the most difficult 
to measure owing to unknown blended lines that prob¬ 
ably cause the quite large point spread, particularly for 
stars with planets. The [Mn/Fe] ratio generally tends to 
increase with the metallicity, differently from other iron- 
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Fig. 8. [X/Fe] vs. [Fe/H] trends using binned average values. The bins are 0.2 dex wide, centred on [Fe/H] 
= —0.6,—0.4,—0.2,0,0.2 and 0.4 for both samples. Triangles are planet-hosts while rhombi are stars without a 
planetary-mass companion. The error bars represent the standard deviation about the mean value. 
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peak elements (see Fig. [5J. In Fig. 0we observe a clear 
change of slope around — 0.2 < [Fe/H] < 0. 

Despite a certain scatter, we note the good agreement 
between our results and those of BOD and SD02. 

Nissen et al. (2000) analysed Mn abundances, in disc 
stars (—1.4 < [Fe/H] < 0.1). These results show an evident 
increase in manganese abundances for the entire metallic- 
ity range considered. Another study (FG98) also exhibits 
a slight linear [Mn/Fe] dependence on iron in metal-rich 
stars. 


6.2.2. Vanadium 


On the subject of vanadium abundances, Fig. 0 (bottom 
right) clearly shows how removing cooler stars (T e gf < 
5000 K) from the data considerably reduces the disper¬ 
sion of points. Although vanadium belongs to the iron- 
peak group, we note that [V/Fe] ratio behaves like an a 
element. Figure |E1 shows that [V/Fe] values in stars with 
planets are systematically ~ 0.10 dex lower than in the 
comparison sample stars. Since planet-hosts are, on aver¬ 
age, hotter than comparison sample stars (see Fig. n and 
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Fig. 10. [X/Fe] vs. [Fe/H] (in dex). Results by Beirao et al. (2005) with the addition of abundances from new 
comparison sample spectra. On the left are the plots for Na, Mg and Al, and on the right the same plots only for 
T e ff > 5000 K objects. Triangles and asterisks represent stars with planetary-mass companions and “single” field stars, 
respectively. The intersection of the dotted lines indicates the solar value. 
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[V/Fe] shows a negative slope with T e g , this effect may 
contribute to the observed difference. 

The same considerable scatter is observed in the stud¬ 
ies by BOD and SD02, but good agreement among trends 
is found in both cases. Contrary to our results and those 
of BOD, SD02 emphasized that the [V/Fe] values of stars 
with planets are about 0.15 dex higher than comparison 
sample vanadium abundances. This mismatch could sug¬ 
gest that vanadium analysis is strongly influenced by an 
NLTE effect, as discussed in Section rm 

There are not many studies regarding vanadium abun¬ 
dances in metal-rich stars without planets. The [V/Fe] ra¬ 
tio remains around the solar value in the plots of FG98 
and CHOO, for [Fe/H] > 0. 

6.2.3. Chromium 

Chromium abundances in our target are constant around 
[Cr/Fe] ~ —0.05 dex (see Fig.0. We note very little scat¬ 
ter for this element (Fig. 0. 

Our observed trend is similar to those already pub¬ 
lished by BOD and SD02: the point distribution is uni¬ 
form, constant and with little scatter. However, the BOD 


data show about 0.05 dex systematically lower chromium 
abundances compared to our results. 

Good agreement is found with the FG98, CHOO, BEN 
and FUL plots. All these results exhibit a [Cr/Fe] ratio of 
around 0. 

6.2.4. Cobalt 

The [Co/Fe] ratio first decreases to the solar value and 
then slowly rises for metal-rich stars (see also Fig. 0. The 
figures also show that the abundances of metal-rich stars 
with planets seems to be lower (~ 0.05 dex) than [Co/Fe] 
values for stars without planets. 

Previous studies of BOD and SD02, have obtained sim¬ 
ilar results for cobalt. In this last paper authors suggested 
that planet-host stars exhibit a slight Co overabundance 
with respect to the solar value (~ 0.15 dex) when com¬ 
pared to stars with no companions. 

Our results resemble those of FG98, while with respect 
to AL04 we obtained a change of slope of around [Fe/H] ~ 
0 . 
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Fig. 11. The same as Fig.|Hlbut for Na, Mg and Al. 

6.2.5. Nickel 

Similarly to calcium and vanadium, the nickel trend shows 
a possible plateau in the range —0.2 < [Fe/H] < 0 and 
hence a slight slope change for metal-rich stars (see Fig. EJ. 

Our trends resemble those observed by BOD, SD02 
and Fischer & Valenti (2005). The [Ni/Fe] ratio seems to 
increase slightly in metal-rich stars. 

In most of the studies considered (see FG98, CH00, 
FUL, EAGLNT) nickel abundances exhibit a uniform and 
approximately constant trend around the solar value. In 
the results proposed by BEN, nickel abundances remain 
constant until [Fe/H] ~ 0.2, where the values are greater 
than the rest of the set of data. A different plot is rep¬ 
resented by AL04: the [Ni/Fe] ratio increases for higher 
metallicity values. 

6.3. Na, Mg and Al 

In Figs ll OlTTI we present the plots for Na, Mg and Al. As 
discussed in detail in the paper by Beirao et al. (2005), 
the [Na/Fe] ratio slowly decreases as a function of [Fe/H] 
until the solar metallicity, then we can observe a change 
of the slope. This behaviour is similar to that shown by 
some refractory abundance distribution (e.g. Si, Ti, Sc). 
Despite the dispersion of the points, the [Na/Fe] values 
are on average below solar for [Fe/H] ~ 0.0. 


The Mg and Al abundances also resemble those derived 
for Si, Ti, Sc. We note that the decrease in the [Mg/Fe] 
and [Al/Fe] values with increasing metallicity in the range 
—0.70 < [Fe/H] < 0 is stronger than for Na. For metal-rich 
stars the abundance distributions for the three elements 
stay at approximately solar value except for a light upturn 
in [Na/Fe] and [Al/Fe] values. 

These results agree globally with previous studies of 
stars with planets (see SD02, GZ00 and Fischer & Valenti 
2005 for Na plots) and we also confirm that Mg and Al 
may be slightly enhanced in the planet-host HD 168746, as 
noticed in SD02. The [Na/Fe] upturn above solar metal¬ 
licity is also found in GZ00 and both EAGLNT, FG98, 
BEN, but not in CH00. With respect to Al abundances, 
other authors (EAGLNT, FG98, BEN) have observed sim¬ 
ilar trends. In the plots proposed by CH00 we note a clear 
upturn at [Fe/H] ~ —0.2,—0.3, as traced by comparison 
stars in our sample. 

Possible differences between the abundance trends of 
the two groups of stars have been observed here only for 
metal-rich stars. It is particularly interesting to note that 
comparison stars continue to show a decrease in [Mg/Fe] 
with increasing [Fe/H], while maintaining the slope ob¬ 
served for lower metallicity, while planet-hosts change the 
slope. An opposite effect is observed for [Al/Fe] where 
the planet hosts seem to have less Al than their “sin¬ 
gle” conterparts. The dependence of [Mg/Fe] on T e g is 
not responsible for this effect (see figs 1 and 2 in Beirao 
et al. 2005) and we cannot rule out that this difference 
is real. Figure QJ]| also shows that removing cooler stars 
(T e ff < 5000A') from the dataset does not change these 
results. GZ00 also suggested a difference in [Mg/Fe] val¬ 
ues between planet-host and single stars. 

7. Discussion and conclusion 

We have determined abundances for nine refractory el¬ 
ements (other than iron) in a large sample of 101 stars 
with planets and in a homogeneous comparison sample of 
94 stars with no known planets. We have also presented 
Na, Mg, Al abundances in 53 new comparison sample 
stars to extend the previous work by Beirao et al. (2005). 
For each element a uniform and independent study of the 
two samples was carried out using atmospheric parameters 
derived from a detailed spectroscopic analysis by Santos 
et al. (2004a, 2005). Abundance ratios [X/H] vs. [Fe/H] 
and [X/Fe] vs. [Fe/H] have been plotted to compare the 
two samples and to try to find differences eventually con¬ 
nected to the presence of giant planets. This study was 
also intended to provide a complete comparison in the 
high metallicity domain, where studies had lacked “sin¬ 
gle” stars with [Fe/H]> 0.1. Furthermore, the data could 
provide clues clarifying the chemical evolution of plane¬ 
tary systems. 

In our analysis we stressed a certain diversity of trends 
for elements of a common origin. On one hand not all the 
refractory elements studied here show the same behaviour; 
on the other hand, the abundance trends of elements com- 


















14 


G. Gilli ct al.: Abundances of refractory elements 


ing from the same nucleosynthesis source are not always 
alike, in contrast with that we expected. 

Our concluding remarks are as follows: 

— Again we confirm that the excess of metallicity ob¬ 
served for planet hosts is not unique to iron. 

— The abundance trends of stars with planets are very 
similar to those traced by comparison sample stars. 
This feature could favour the primordial hypothesis 
to explain the metallicity excess in stars harbouring 
planets. In any case, some elements (e.g. Mg, Al, V, 
Co) show certain differences in the behaviour of abun¬ 
dances of stars with planets and “single” stars, in the 
higher metallicity range. We thus do not exclude the 
possibility that the presence of a planet might influence 
the composition of certain elements in the atmosphere 
of metal rich stars. 

— Good agreement was found with both previous pub¬ 
lished studies on abundances in stars with extrasolar 
planets and most studies on metal-rich stars of the 
Galactic disc. One might suggest that the observed 
trends are simply a conseguence of Galactic chemical 
evolution, with no particular mechanism linked to the 
presence of a planet. To this end, only the calcium 
abundances show a different trend when comparing 
studies on stars with planets (BOD, GZOO and this 
article) and most chemical analysis on Galactic disc 
stars (EAGLNT, FG98, AL04, FUL). 

In the future gathering new homogeneous abundance 
values of other elements with a wide range of conden¬ 
sation temperatures Tq will be of great importance. For 
example, a detailed comparison of [X/Fe] abundances of 
volatile and refractory elements is currently in progress . 
This study will give us the chance to discuss the relative 
importance of differential accretion (e.g. Gonzalez 1997; 
Smith et al. 2001; Sadakane et al. 2003) in stars harbour¬ 
ing extrasolar planets. 
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